Cosmic string evolution with a conserved charge 
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Cosmic strings with degrees of freedom beyond the standard Abrikosov-Nielsen-Olesen or Nambu- 
Goto strings are ubiquitous in field theory as well as in models with extra dimensions, such as string 
theoretic brane inflation scenarios. Here we carry out an analytic study of a simplified version of 
one such cosmic string model. Specifically, we extend the velocity-dependent one-scale (VOS) string 
evolution model to the case where there is a conserved microscopic charge on the string worldsheet. 
We find that whether the standard scale-invariant evolution of the network is preserved or destroyed 
due to the presence of the charge will crucially depend on the amount of damping and energy losses 
experienced by the network. This suggests, among other things, that results derived in Minkowski 
space (field theory) simulations may not extend to the case of an expanding universe. 
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I. MOTIVATION 

Cosmic strings 0, , line-like topological defects that 
can be produced in cosmological phase transitions, pro- 
vide a valuable tool for constraining models of the early 
universe. Having a wide range of potential observational 
signals Q, which depend directly on their microphysi- 
cal properties, they can be used to constrain high-energy 
physics parameters from cosmological observations. In 
particular, the magnitude of their observational signa- 
tures is mainly determined by the energy scale of the 
corresponding symmetry breaking transition. Interest 
in the study of cosmic strings has been renewed in re- 
cent years [|| with the realisation that string-like de- 
fects are generic in a wide range of fundamental mod- 
els, including Supcrsymmctric Grand Unified Theories 
Q and Brane Inflation models in String Theory 
In the latter case, there is even the exciting possibility 
of constraining (within a given model) fundamental pa- 
rameters like the string coupling, string scale, and warp- 
ing /compactification scales from cosmological observa- 
tions [itniini. 

In order to be able to make this quantitative link be- 
tween late time observations and early universe physics, 
one must therefore be able to model the evolution of 
string networks from the time of their formation until 
they are observed. Work in this direction has mostly 
focused on the simplest, structureless strings of the 
Abrikosov-Nielsen-Olesen type and their corresponding 
Nambu-Goto thin string approximations. In recent years, 
more complex models have been developed (both on the 
numerical [l2l - fl5| and analytical [l6l - [Tc| fronts) , allowing 
for the formation of zipping interactions between strings 
of different tensions, in order to describe the richer nature 
of the so called FD-nctworks that appear in Brane 
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Inflation models. 

However, even networks of a single type of string can 
have non-trivial internal structure, generally carrying ad- 
ditional degrees of freedom on the string worldsheet. This 
is the generic situation in models with extra dimensions, 
where there is a proliferation of scalar fields that can cou- 
ple to (and condense on) the strings. Cosmic strings with 
additional worldsheet degrees of freedom (scalar charges, 
currents, fermionic zero-modes) have been extensively 
studied in field theory and supergravity, see for example 
[20U25I ]. One important result that emerged from these 
studies was the realisation that the presence of currents 
can lead to the formation of stable cosmic string loops, 
known as vortons [22l [26j , which could therefore domi- 
nate the energy density of the universe. This has allowed 
to obtain strong bounds on models producing vortons 

To date, the cosmological evolution of such string net- 
works remains largely unexplored. It is therefore desir- 
able to understand how these additional degrees of free- 
dom can be described macroscopically, and how their 
presence affects the behaviour and cosmological conse- 
quences of the corresponding string networks. In a recent 
paper [3(| , the evolution of semilocal strings [U was de- 
scribed analytically and it was found that the spectrum of 
scaling solutions has a much richer structure than in ordi- 
nary cosmic strings. The presence of worldsheet degrees 
of freedom would also modify significantly the structure 
of scaling network solutions and can be expected to ob- 
struct scaling, possibly leading to string frustration. 

Here, we focus on the case of cosmic strings with a con- 
served charge li ving on the string worldsheet. We extend 
the VOS model |32H35j to describe this case analytically, 
and study the effect of the charge on the evolution of the 
string network. A recent related example, for the case of 
domain walls, comes from the work of Battye et al. [36[ 
(but see also [37| for some caveats), who suggest that 
with a conserved, homogeneous and localised Noether 
charge the number of domain walls does not scale in the 
usual way, in which case they could provide some con- 
tribution to the dark energy suggested by cosmological 
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observations. These results are based on relatively small 
numerical simulations in 2 + 1 dimensions and without 
cosmological expansion, so extrapolating them to a cos- 
mological context requires some care. Nevertheless, they 
highlight the need for a better macroscopic description of 
these processes. As we shall see, our analysis will provide 
some possible explanations for these results. 



in Ref . [28J . A comprehensive description of the general 
formalism can be found in 13811 . 



Generic chiral models 



II. MICROSCOPIC MODEL 

Here we review the basics of the microscopic model 
that we will be using. This model is described in detail 



For strings carrying conserved currents, we consider 
the Witten-Carter-Peter (chiral) model, which implicitly 
makes use of the fact that in two dimensions a conserved 
current can be written as the derivative of a scalar field. 
The dynamics is described by the Witten action 
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where the four terms are respectively the usual Nambu- 
Goto term (//q being the string tension and 7 the pullback 
of the background metric g on the worldsheet), the inertia 
of the charge carriers described by the scalar </>, the cur- 
rent coupling to the electromagnetic potential A^, and 
the kinetic term for the electromagnetic held. World- 
sheet indices are denoted by a,b <E {0,1} and we will 
take <t° to be the timelikc coordinate, while a 1 = a will 
be spacelikc; e ah is the alternating tensor in two dimen- 
sions. Note that this action applies to both the bosonic 
and the fermionic case. 

We are essentially interested in the chiral limit of this 
model, that is (taking dot/prime to denote differentiation 
with respect to the timelike/spacelike worldsheet coordi- 
nate cr /a): 



and choose the standard gauge 



cr° = T , X • x' = . 



in which the scalar e becomes: 



(5) 



2 \ 1/2 



1-x^ 
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Then, in the chiral limit, introducing the simplifying 
function $ defined as 



b 12 = e V . 



where e is the scalar 



(2) 



(3) 



giving the string energy per unit coordinate length. Let 
us now consider an FRW background 



$(0) 



(7) 



/z a 2 (l - x 2 ) ' 
the microscopic equations of motion take the form [28j : 



[e (1 + $)]•+ f x 2 = $'-2-e$. 

id a 



(8) 



ds 2 = a 2 (dr 2 - dx 2 ) 



(4) 



and 



e(l + $)x+— (l-x 2 )x = 

id 



(1-$) 



$ + 2-$ I x' + 2$x . 

a 



(9) 



r 



where for simplicity we have introduced the damping describing Hubble friction (H = a 1 (da/dt) = a/a 2 ), but 
length also allowing for an additional friction mechanism with 



= a I 2H + -r 



(10) 
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characteristic lengthscale if. 
Two generic useful relations are 



and 



1 -x 2 



1 -x 2 



(11) 



(12) 



The curvature vector cPx/ds 2 , with ds = V x' 2 dcr, satis- 
fies: 



fx _ _ (x' ■ x")a 
ds 2 x' 2 



r'4 



(13) 



Thus, a curvature radius, i?, can be defined locally via 



e(l-x 2 ) 



X • X X ■ X 



r>2 



r>2 



R 



(x-u), (14) 



where we have introduced a unit vector u in the direction 
of the curvature vector. 

The worldshcet charge and current densities are re- 
spectively given by 



qe<p, 



and 



(15) 



(16) 



while the total energy of a piece of string is given by 

E = fi a J (1 + $) eda = E s + . (17) 

We can immediately interpret the energy (|17p as being 
split in an obvious way into a string component and a 
charge component. Defining a macroscopic charge as the 
average of $ over the string worldsheet, we then have: 



= (*> 



_ / §eda E<s> 



J eda 



(18) 



This interpretation will be relevant below. 

Introducing the network string density p s such that 
E s oc p s a 3 and defining the correlation length £ by 



_ Mo 

Ps - £ 2 , 

the evolution equations have the form 



E ± = P 1 + 3 a 
E s p s a 



£ da 
-2^ + 3- = - 
£ a a 



and 



V = (x 2 ) , iw = (x • x) 



(*>• 



(19) 

(20) 
(21) 

(22) 



Two further (not independent) useful relations are 

E<s> _ Q E s 
E$ Q E s 



(23) 



and 



E E° Q pa La , 

- = - i + --^ 7 r = -+3- = -2- + 3-, 24 
E E s 1 + Q p a L a y ' 

where p oc Ea~ 3 is the energy density and the lengthscale 
L is defined by: 

Mo 



L 2 



(25) 



With these, one can in principle proceed to average this 
model. An additional difficulty which is absent in the 
case of Nambu-Goto strings is the appearance of a term 
proportional to (x' • x"). This factor is not expected to 
be zero even though (x' • u) is, see equation ([TB"]). This 
will be discussed in more detail below. 

Finally, note that in the case of the VOS model for 
plain Nambu-Goto strings one assumes that the network 
has a single characteristic length scale, so that R = L = 
£. This is no longer true in the charged case, but we will 
still assume that R = £, while L is now only a measure 
of the total energy in the network. 



B. Conserved microscopic charge case 

In our particular case we have 

Pw = oecj) = q(j) = const. 

and therefore 



$ = 



Vo 



(26) 



(27) 



with tpo being a constant. Now, by simple differentiation 
one finds that $ evolves as follows 

a x • x" 

$ + 2-$ = 2$ — (28) 

a x' A 



and 



$' + 2$- 



r'2 



0. 



(29) 



The conserved charge assumption simplifies some of 
the above equations. If one expands Eq. (jlip by succes- 
sively substituting in Eqs. pfe|) and then Eqs. (|28H29|) 
one obtains 



(30) 
(31) 



and inserting this in Eq. (fT2"j) yields 
e' = 0. 



An alternative way to see this is to substitute Eq. © 
into Eq. (JT2J) and then use Eqs. (|28«29|) . This leads to 

thus if $ 7^ we must have e' = 0. 
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III. INTERLUDE: NAMBU-GOTO STRINGS IN 
EXTRA DIMENSIONS 

The above formalism describing strings with currents 
is also analogous to the effective description of struc- 
tureless, Nambu-Goto strings evolving in the presence of 
compact extra dimensions. Indeed, integrating out the 
extra dimensions gives rise to worldsheet scalars which 
geometrically correspond to string positions in the inter- 
nal manifold and behave like currents in the low energy 
worldsheet theory. 

The evolution and dynamics of strings moving in space- 
times with extra compact dimensions have been studied 
in Ref. j^. Consider the following 'augmented' FRW 
metric 



ds 2 = dt 2 - a{t) 2 dx 2 - b(t) 2 d\ 2 , 



(33) 



where the internal manifold (coordinates 1) has been 
taken for simplicity to be toroidal with scalcfactor b(t). 
Now let the vector l(r, a) denote the string position 
in this internal manifold. In the simplest static case, 
b(t) = 1, the string energy is: 




where we have expanded to linear order in l' 2 /(a 2 x' 2 ) 
and i 2 /[a 2 (l-x 2 )]. 

The analogy to the above discussion is apparent: the 
first factor in the last line of (fM)) is simply e, while the 
terms in the second factor are of the form of equations 
and (HJ. Indeed, setting i 2 /[a 2 (l - x 2 )] = $, and 
considering the analogue of the 'chiral limit' 4>' 2 = e 2 c/> 2 
discussed above, we then have in this case 



The embedding fields 1 = l m (m = 1,...,D), describing 
the string positions in the D compact extra dimensions, 
are worldsheet scalars which in general have richer dy- 
namics than in action (fTJ). In particular, their kinetic 
structure is governed by the relativistic action Nambu- 
Goto action in 4 + D spacetime dimensions: 



S = - 



Mo / ^-dct (gffdaxrdbx'+g^dami™) , 



(38) 

where g^J is the (4-dimensional) spacetime metric and 
ifp the metric in the D compact dimensions. Factoring 



out 7^ = g^daX^dbX^ yields: 



Mo / \/-det 7 $ y'det^ + 7 - dj'hl™ g£l) . 



S 

J (39) 
that is, the general action for the worldsheet scalars has 
non-trivial kinetic structure. However, for small I and /', 
the second factor in (|51))) can be linearised, and with the 
additional 'chiral' condition (|35[) one recovers the first 
two terms of the Witten action (JXJ) , for D scalars, l m , 
subject to (|35j) . 



IV. LOOP SOLUTIONS 

Simple exact loop solutions can be used to further our 
understanding of the role of the various terms on the 
evolution of the strings, and also determine the accuracy 
of averaged (macroscopic) quantities in relation to micro- 
scopic quantities. The simplest loop solution is a circular 
one, which has the form 



x = r (r) (cos 6, sin 6, 0) 



(40) 



with 



12 



r'2 



(1-X 2 ) 

and equation (fM)) becomes 



(35) 



and 



VT^T 2 



(41) 



e e ~e(l + $) 

in direct analogy to equation (|17l) . In addition, the equa- 
tion of motion for e e 



^ = -^[l + x 2 -x' 2 /e 2 ] 



(36) 



with e e ~ e(l + $) becomes: 



[e(l + $)]• = 



-2-e(x 2 + $) + 0($ 2 ), (37) 
a 



which is equation © for $' = to linear order in <f>. 

Note that this analogy only becomes quantitative in 
the limit where l' 2 /(a 2 x' 2 ) and l 2 /[a 2 (l — x 2 )] are small. 



Mo a 



eda . 



(42) 



For this ansatz Eqs. (JSJ and © arc equivalent, and 
either of them can then be written 



(1 + *)r + (1 - f 2 ) 



d , ,1 
2-f + 1 - $)- 

a r 



, (43) 



while the evolution equation for $ has the form 



2$ | H + - 

r 



(44) 



which is in fact trivial given the definition of $ in (1411) . 
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We start by examining the case with no expansion 
(H = 0). The evolution equations, Eqs. fl£J) and (0), 
now simplify to 



e(l + <J>) = const. = 1 



and 



r + (1- $ 2 )r = 0. 
Alternatively, they can be written as 



r 2 = l-2^-4^ 



and 



r + r 



4 



(45) 
(46) 

(47) 

(48) 



These have periodic (oscillatory) solutions, but unlike 
the case of circular loops in the plain Nambu-Goto case, 
the presence of a charge ensures that these loops never 
collapse to zero size and that their microscopic velocity 
is always less than unity. Indeed, the maximum velocity 
is 

fi ax = 1-4$. (49) 
The limiting behaviour corresponds to a static solution, 

r = ; (50) 



1 



r = ip = 



in this solution the energy of the loop is equally divided, 
half of it being in the string itself and the other half in 
the current. In the general case, the average loop velocity 
(over one oscillation period) is 



1 



(r^) = -(l-4^), 



(51) 



and the average fraction of the loop's total energy in the 
current is 



~E~ 



<Po ; 



(52) 



naturally the rest of the energy corresponds to the (bare) 
string. Finally the macroscopic charge of these loops, 
given by Eq. (fl"8|) . is 



Q 



1 



¥>0 



(53) 



notice that this is unity in the static limit (<po = 1/2). 
The static solution we discussed above still exists in an 
expanding universe. In this case it corresponds to 



ar = ip , r = 0; 
such a loop will have a physical radius 

n.r 

Rphys 



(54) 
(55) 




FIG. 1: The time evolution of a particular loop solution in 
an expanding universe. The solid (blue) curve shows r, while 
the dashed (red) one shows the force (po/r 3 due to the charge, 
Eq. (|48[1 . As the loop shrinks and r approaches zero, the 
force provides a 'boost' causing it to expand again. 



Solving the differential equation for r and integrating 
over one period, we can obtain 



(x 2 ) 2 
<x 4 ) 



(56) 



which is the same as in the case without charge. This 
result is important because the fact that it is always of 
order unity allows us to, every time a microscopic quan- 
tity such as x 4 appears in an expression inside an integral, 
average it to v 4 to a good approximation. In the general 
case with expansion, Eq. (|43[) cannot be solved analyti- 
cally, but numerically we found that this result still holds 
to a good approximation for a wide range of initial condi- 
tions. Fig.fTJshows the behaviour of a particular solution, 
for which we find 



0.6622. 



(57) 



As can be seen in Fig. [TJ the charge in the loops keeps 
them from reaching r = by providing a 'boost' when r 
is approaching zero. It then causes r to grow again (to a 
lower value than before), while it returns to its previous 
value. Naturally, these simple loop solutions satisfy Eqs. 
(|30H3ip . and indeed one can check that this is still the 
case for Kibble- Turok [4(| and Burden [4l[ loops. 



V. MACROSCOPIC EQUATIONS 

We can now proceed and look at the averaged evolution 
equations in our conserved microscopic charge case. 



A. General dynamical equations 

The total energy of the string is given by Eq. (fT7]) and 
therefore we can define two characteristic lengths for the 
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string, as in section Hi At the usual correlation length £ 
associated with the string energy 



E s = fia / eda oc p s a 3 



(58) 



through Eqn. ([T9")) , and the lengthscalc L associated with 
the total energy E oc pa 3 , sec Eqs. (|T7|) . (|25|) . Taking 
the time derivative of the previous equation, one easily 
finds the evolution equation for £ 



2H 1 



2Q kv 



l + Q J 1 + QR \e(l + $) 



(59) 



Here, k is the so-called 'momentum parameter' quantify- 
ing the average angle between the curvature vector and 
the velocity of string segments in the network. It thus 
provides a measure of the small-scale-structure on strings 
and can range from (wiggly strings in flat space) to or- 
der unity (smooth strings) [34j . 

Correspondingly, the evolution equation for L is 



L 
> 

' L 



2H 1 







l + Q 



e(l + $) 



(60) 



As was defined above, Eqn. (fT5|) , the macroscopic charge 
Q is the averaged microscopic charge (<&), which is also 
the ratio between the string energy E s and the 'charge' 
energy E$. Differentiating (fT8)| . we find that its evolution 
equation is 



Q 



= 2 



kv 
~R 



- H 



Finally, the string velocity is defined by 
2 2 J x 2 eda 



(61) 



(62) 



and taking time derivatives on both sides, we arrive at 
the evolution equation 



1 



1 + 



A(l_Q)_ 2Ft , + i± 

R v 



e(l + $) 



(63) 

However, equations ([59^1 . (J5DJ) and (|5T|) are not indepen- 
dent: there is a consistency relation between £ and L 



E s 



E 



i = L^Jl + Q 



1 + 

Therefore, the equations are related by 



2 ^ = 2 i-- 

L £ 1 



(64) 



(65) 



which is verified, so the three equations are consistent. 



B. Ansatz for the charge gradient 

As previously mentioned, in order to proceed we now 
need to deal with the (x' • x") term coming from 
Referring to Eq. (|29[) . dimensional analysis suggests an 
ansatz of the form 



e(l + *) 



v 2Q 
! R1 + Q 



(66) 



where s is (at least, to a first approximation) a constant. 

Using (|66[) and noting our earlier identification R — £, 
our evolution equations become: 



2 Hn i+ — 



20 (fc + s)t> 
"+Q 



2Q 



i _ / v 2 + Q N 

2 I- 2ff + TToV ' (1 + O) 3 / 2 l 



Q_ fkv 



if 



i_^ 

l + Q 



(1 - Q(l + 2s/fe)) - 



(67) 



(68) 



(69) 



(70) 



We will assume a critical-density universe with generic 
expansion rates of the form 



a (x r 



H 



A 

7 



and look for scaling solutions of the form 

£ = & a 

(or an analogous law for L), 

v = vot 13 , 

and 

= Ooi 7 ■ 



(71) 

(72) 
(73) 
(74) 



Note that causality implies a < 1 and the finite speed of 
light implies ft < 0. Furthermore, our discussion of loop 
solutions shows that v — > 1 is not a physically allowed 
solution for these networks. 



VI. NO CHARGE LOSSES 

In this section we assume that there are no macro- 
scopic charge losses. (The case with charge losses will be 
discussed in the following section.) We will separately 
consider the cases with and without energy losses due to 
loop production. 
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Whether or not we have loop production, the evolu- 
tion equation for the macroscopic charge Q is given by 
Eq. (|69)) . and we can start by studying this. There is 
a trivial but unphysical (refer to discussion after Eqn. 
(|59[1 ) solution if k = 0, with £ oc L oc a and v oc Q oc a~ 2 , 
which can therefore be ignored. In the realistic case k ^ 
there can in principle be two kinds of solutions: 

• Decaying charge solutions, with 



• For A = 2/3, corresponding to the matter- 
dominated era, 



= 0, 



a = 1 , Co = -A^o 



-o 2 



1 + | 



(81) 



(82) 
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-2A, /3<a-l 



(75) 



for these solutions not only does the charge decay 
(as Q oc a" 2 ) but velocity will necessarily decay as 
well. 



standard solutions with 



P = a — 1 , 7 



-2A + 2A; 



vo 
Co' 



(76) 



here we have used the term 'standard' referring to 
the fact that linear scaling solution (with a = 1 
and (3 = 0) is of this form, although a priori there 
is no guarantee that this solution will exist with a 
constant (non-zero) charge. Also note that in this 
branch of solutions we may at least in principle have 
growing, constant, or decaying Q. 

We can now study the entire system of equations in 
the cases with and without energy losses. 



and 



7 = 



Ps 



p 



Pcrit i + Qo 



16tt 
3fc2 



1 - 



2s 
~~k 



(83) 



G/i ; 
(84) 

in this case the macroscopic charge is still a con- 
stant, but the additional dilution caused by the 
faster expansion rate is enough to ensure that the 
energy density of the network is a constant fraction 
of the background one. In other words, for this par- 
ticular expansion rate we have a generalised linear 
scaling solution, with £ (and L) growing as fast 
as allowed by causality, in which the RMS veloc- 
ity and the macroscopic charge are constant (and 
larger charges leading to smaller velocities, as was 
seen in the loop solutions). 

For fast expansion rates, A > 2/3, 



A. Without energy losses 

In this case we obtain the following three scaling rela- 
tions: 



For slow expansion rates, A < 2/3, 

3 . 1 t kv 
a=-X<l, £o„ 



/3 = a- 1< 



2s 

{o = I 1 + — 



and 



oc oc i 2 ~ 3A 



Pcrit Pcrit 



(77) 



(78) 



(79) 



(80) 



here we have a constant charge, which gradually 
slows the strings (making v —¥ 0), although the 
evolution is still faster than conformal stretching 
(which corresponds to a = A). As a consequence, 
both the energy density in the strings and the total 
energy density in the network grow relative to that 
of the cosmological background. 



k 2 



(3 = 0. 



4A(1 - A) 



1 - A 



4- 6A < 



(85) 



(86) 



(87) 



and 



p 



32?r 1 - A 



Pcrit Pcrit 3fc 2 A 



Gp ; 



here the additional Hubble damping makes the 
macroscopic charge decay, and we therefore end up 
with the solution for ordinary Nambu-Goto strings 

[II, 11, IH • 

Therefore, which of these solutions is the attractor for 
the network's evolution depends on the universe's expan- 
sion rate, which is parametrized by A. (In principle, one 
also finds a fourth solution which has vq = 1, but, ac- 
cording to our discussion in the loops section, this is un- 
physical for charged strings.) By numerically solving the 
evolution equations one can confirm that the above solu- 
tions are indeed the attractors in the relevant parameter 
ranges; Fig. [5] shows examples for the two scaling regimes. 
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FIG. 2: Illustrating the two scaling solutions in the case with- 
out energy losses. The plots show the relevant macroscopic 
quantities divided by the time dependencies inferred in this 
section's analysis, so that asymptotically constant values im- 
ply convergence to the expected scaling solution. Top panel: 
Case A = 1/2; from top to bottom the lines depict £/t 3//4 , Q 
and vi 1 / 4 . Bottom panel: Case A = 3/4; from top to bot- 
tom the lines depict £/t, v and Qt 1//2 . In both cases we have 
assumed k = 0.8 and s = 0.01. 



For slow expansion rates, A < 2/(3 + c/fc), 



P = a- 1 < 



7 = 0, Qo = 1 



2s 
~k 



and 



Pa 



OC 



Pcrit Pcrit 



— oc t 2 - {3+c/k)x ; 



(91) 



(92) 



(93) 



(94) 



loop production is an additional energy loss mech- 
anism, and therefore the correlation length now 
grows faster than in the c = case. Similarly, 
string velocities decrease more slowly, and the net- 
work's density grows more slowly relative to the 
background one. 



• For an intermediate expansion rate, A 
c/fc), 

a = 1 , £ = iu (3fc + c) 



= 0, 



-o 2 



1 + l 



7 = 



2/(3 + 
(95) 

(96) 
(97) 



and 



B. With energy losses 

The loss of energy through loop production is described 
by a new term in the evolution equation for £ (and, con- 
sequently, one for L), which, using standard arguments 
[33, K33l [42f , one can write as 



(89) 



This introduces the parameter c quantifying the efficiency 
of producing loops. Similarly, using Eq. (|65[) and assum- 
ing there is no charge loss, we find the equivalent term 
for L, 



4 



i 



(90) 



The analysis of scaling solutions can now be repeated, 
and one finds solutions that generalise the above ones by 
including an additional dependency on the loop chopping 
efficiency c: 



P 



Pcrit 



!0 Pcrit 



167T 
3fc2 



+ ¥j Gm; 

(98) 

now the expansion rate for which this solution oc- 
curs decreases, with the smaller expansion rate be- 
ing compensated, for the purposes of the dilution 
of the network's energy density, by the process of 
loop production. (For c = this solution exists for 
the matter era, and c = fc will make it occur in 
the radiation era.) Interestingly, the ratio of the 
string and background energies is exactly the same 
as before — in other words, it does not depend on 
the value of c. 



• For fast expansion rates, A > 2/(3 + c/fc), 

fc(fc + c) 



/3 = : 



4A(1 - A) 



1 - A fc 
A fc + c 



(99) 



(100) 
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FIG. 3: Illustrating the two scaling solutions in the case with 
energy losses. The plots show the relevant macroscopic quan- 
tities divided by the time dependencies inferred in this sec- 
tion's analysis, so that asymptotically constant values imply 
convergence to the expected scaling solution. Top panel: 
Case A = 1/2; from top to bottom the lines depict £/t 0,822 "', 
Q and vt 0,178 "'. Bottom panel: Case A = 3/4; from top to 



bottom the lines depict £/t, v and Qt 
have assumed k — 0.8, c = 0.23 and s — 0.01. 



In both cases we 



and 



k--(3k + c) 



32tt 1 - A 



Peril Peril 3fc(fc + c) A 



< 



(101) 



(102) 



this is exactly the VOS linear scaling solution for 
Nambu-Goto strings (32TJ35l | . with an added predic- 
tion of a particular decay law for the charge (which 
depends on the cosmological expansion rate) . 

Fig. [3] illustrates the two branches of the solution. It 
can be trivially seen that in the limit of no loop produc- 
tion, c — > 0, these solutions are the same as in the case 
without energy losses, Eqns. (|77|) - ([55|) . As before, there 
is a fourth solution with v = 1 which is not allowed for 
charged strings. The analysis so far also highlights the 
point that the parameter c is important for determining 
whether charge survives on the strings or decays. This 
motivates a discussion of possible charge losses, which 
will be done in the next section. 



VII. MACROSCOPIC CHARGE LOSSES 

In the previous section, in order to derive the energy 
loss term in the evolution equation for L, it was assumed 
that there was no charge loss. Now this assumption will 
be dropped, and the resulting scaling relations analysed. 
This section is somewhat more phcnomenological than 
the rest of the article, since it will rely on simplifying 
assumptions for how the charge may be lost, but our goal 
is simply to develop an intuitive picture for the possible 
role of charge loss mechanisms on the evolution of the 
network. 

We will again assume an energy loss term of the form 



r 2 t 



(103) 



but this time let us say that L gets a different (in general) 
loss term 



E 
E 



-fc 



L 



L 
> 

' L 



fc 



L ' 



(104) 



where / is an arbitrary function, possibly of the velocity, 
charge and correlation length. We are therefore assuming 
that any such charge losses are related to the network's 
intcrcommuting and loop production mechanisms. 

In the previous section, we assumed there was no 
charge loss, implicitly using 



i 



(105) 



Instead, we will now leave / free and obtain the evolution 
equation for Q from the above assumptions together with 
Eq. USSJl. The result is 



Q = 2Q(— - H) + -(1 + Q) [1 - fy/T+Q) , (106) 

or an analogous equation in terms of L. With the choice 
(|105p for / we trivially recover the results of the previous 
sections. We can now look for scaling relations as before, 
checking how the results depend on the choice of /. 

We start by noting that, in the standard VOS-model 
case (without a charge), we would not expect charge to 
be created. Therefore, Q = when Q = 0, implying 



2y(W) = 0. 



(107) 



It then follows that, if the function / is constant, it has 
to be equal to unity. Otherwise it must be dependent on 
the charge (and possibly other quantities), and become 
unity when the charge drops to zero, 



f = f(Q,v,0 



(108) 



as was the case for the function (|105p we used in the case 
without charge loss. 

Let us therefore consider the / = 1 case. The gener- 
alised scaling laws now become: 



10 



• For slow expansion rates, A < (2k — Wc) / [3k + (1 

W)c], 



3 

a= 2 



1 + (1 - W)c/3k 



1 - Wc/2k 

P=a-1<0 



A<1, ( = j(k-Wc/2) 



1 = 0, Q 



2s 

1 + T 



and naturally 

Ps 



OC 



P- oc t 2 ^ 



Pcrit Pcrit 



(109) 
(110) 

(111) 
(112) 



(where for simplicity we have kept a in the last 
expression); we have also introduced 



W = (l + Qo 1 ) (yi + Qo-i) 



(113) 



which is always positive and behaves as W = 
0.5(1 + Qo) in the limit Q — > and as W oc VQo 
for Qo ~^ 00 ■ Notice that the scaling exponent for 
£ now has an explicit dependence on the charge, 
which was not the case without charge losses. 

As one would expect, increasing the scaling value of 
the charge pushes the value of the maximal expan- 
sion rate where this regime holds to larger values, 
whereas the scaling exponent a decreases. Simi- 
larly, string velocities decrease faster, and the net- 
work's density grows faster relative to the back- 
ground one. 

In principle, as one makes W progressively larger, 
the scaling exponent a becomes closer to A, which 
corresponds to the conformal stretching case. How- 
ever, in practice we do not expect this to occur, 
since it is clear from its definition that W should 
be a small parameter in realistic (cosmological) de- 
fect networks. 



• For an intermediate expansion rate, A 
Wc)/[3k+(l-W)c], 

a = 1 , Co = \vo [3k + (1 - W)c] 



= 0, v 2 = 



1 



2 - Wc/k 



1 - 



2s 



7 = 



and 



Ps 



1 



p 16ir 

Pcrit 1 + Qo Pcrit 3fc 2 



2s 

l - Qo ( l + -r 



(2k - 
(114) 

(115) 
(116) 



Gp. 
(117) 



As in the previous solution there is now an explicit 
dependence on the amount of charge, and the ex- 
pansion rate for which this solution occurs increases 
with the charge. This is to be expected: a larger 
charge makes scaling harder, requiring more energy 
losses (from the damping due to the Hubble ex- 
pansion, or from losses due to loop production) to 
counteract it. In fact the expansion rate for which 
this solution exists would approach A = 1 as the 
scaling value of the charge Qo becomes arbitrarily 
large, although as we already pointed out we do not 
expect this to occur in practice. The ratio of the 
string and background energies is still exactly the 
same as before — in other words, it does not depend 
on the value of c or W. 

• For fast expansion rates, A > (2k — Wqc) / [3k + (1 — 
Wo)c], 



2 k(k + c) 

a = 1 > Co = 4A(l-A) 



/3 = 0. 



W c A 



1 - A k 
A k + c 



[3fc+(l-W )c] 



(118) 



(119) 



and 



Pcrit 



32tt 1 - A 



p cr it 3k(k + c) A 



Gp ; 



< (120) 



(121) 



which is again the VOS linear scaling solution for 
Nambu-Goto strings, now with a faster decay law 
for the charge (which is obvious since we have ex- 
plicit charge losses). For this solution we used the 
notation Wo to indicate the value of W in the limit 
Qo 0, that is W = 1/2; the reason for this choice 
will become clear below. 

Again it is easy to check that if we set c = and/or 
W = we recover the results of the previous section. 
However, note that the transition between the second 
and third solutions will now depend on the amount of 
charge loss. The expansion rate of the second solution 
coincides with the minimum expansion rate of the third 
solution for 



W = W = I 



(122) 



which corresponds to the Qo = limit. 

Finally, it is interesting to discuss what happens in 
the more general case where / ^ 1. In the absence of 
compelling arguments suggesting a particular form for / 
(other than the case without charge losses which we stud- 
ied in the previous section) we will consider a linearised 
form, that is 



f(Q) = l + wQ, 



(123) 
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charge, with the defect velocities decreasing and the net- 
work correlation length growing more slowly than allowed 
by causality. However, even in this case, the network docs 
not generically frustrate: the correlation length evolves 
more slowly than £ oc t but faster than £ oc a (conformal 
stretching); only in the limit of arbitrarily large charge 
would the conformal stretching behaviour be reached. 

The rate of expansion at which the transition from 
charge domination to linear scaling occurs will depend 
on the dynamical processes acting on the network. A 
larger charge makes scaling harder, but cosmological ex- 
pansion provides a damping term which tends to dilute 
the charge and overcome its influence on the dynamics. 
Further energy loss mechanisms facilitate scaling. Specif- 
ically, this critical expansion rate is given, in terms of the 
model parameters identified in the previous sections, by 



FIG. 4: The value of the expansion rate A at which the at- 
tractor for the network's evolution changes from a fixed charge 
solution to the standard linear scaling one, as a function of 
the parameters w and c/k. 



X = 



2 1 - (0.5 + w)c/2fc 

3 l + (0.5-w)c/3fc 



(126) 



for to real with \w\ < 1. The rationale for this is that 
in realistic networks in cosmological contexts the charges 
are likely to correspond to a small fraction of the overall 
energy density of the network. 

Repeating the analysis we find that the above solutions 
still hold, provided that one extends the definition of the 
parameter W to 



W = (1 + Qo 1 ) [(I + wQ )Vl + Qo- 1 
In the limit Qo ~~ * this now behaves as 

w =(l + w ) (i + Qo); 



(124) 



(125) 



we trivially recover the behaviour in the w = (constant) 
case, but we also see that in this limit W will vanish if 
w — —1/2. This is interesting because that choice of w 
corresponds to the linearised version of / = 1/^1 + Q, 
which as we argued is the case of no charge losses. This 
shows that the above analysis is self-consistent. 



VIII. DISCUSSION AND CONCLUSIONS 

We have extended the velocity-dependent one-scale 
string evolution model to the case where there is a con- 
served microscopic charge on the string worldsheet. Wc 
find that there are two possible regimes for the evolution 
of the network. When the expansion rate of the uni- 
verse is fast enough the macroscopic charge will decay 
and the attractor solution is the standard linear scaling 
regime, where the network is losing energy as fast as is al- 
lowed by causality. When the expansion rate is relatively 
slow, the attractor solution has a constant macroscopic 



Without any losses the exponent is 2/3 (corresponding 
to the matter-dominated era) , and its dependence on the 
model parameters is illustrated in Fig. 2) 

While our study was purely analytical, we should note 
that some numerical work already exists which, at least 
qualitatively, is in agreement with our results. Specifi- 
cally, the aforementioned results of (3(| H3] , which sug- 
gest violations of scaling, rely on simulations which are 
nominally done in Minkowski space, although they have 
some effective numerical damping which probably mim- 
ics a fairly small (but non-zero) expansion rate. This 
therefore agrees with our results for low expansion rates, 
if one assumes that there are some effective charge losses. 
In that case our results predict that the scaling exponent 
for £ will decrease as the charge increases. However, we 
would further suggest that in a cosmological setting (that 
is, with a faster expansion rate), the network may still 
scale. 

Interestingly, something like that kind of behaviour 
seems to emerge in more recent work [43[ , in which appar- 
ent scaling deviations in Minkowski space (for a different 
model, again in two dimensions) are in fact eliminated by 
increasing the amount of damping. A more quantitative 
comparison between our analytic solutions and these nu- 
merical results would certainly be instructive but is not 
possible at this stage: it would require knowledge of the 
behaviour of the defect velocities in the simulations, and 
that is not provided in any of these works [36|, H3, EH . 

In any case, our analysis highlights the point that re- 
sults derived in Minkowski space may not necessarily ap- 
ply to the case of an expanding universe; this has already 
been discussed, for example for the case of wiggly cosmic 
strings in [44[. Finally, our results are also relevant for 
the case of superconducting strings and vortons, but this 
case is left for future work. 
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